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The cross section of the process e + e — » tv + tt tv° was measured 
in the Spherical Neutral Detector (SND) experiment at the VEPP-2M 
collider in the energy regiony's below 980 MeV. This measurement was 
based on about 1.2 x 10 6 selected events. The obtained cross section was 
^ , analyzed together with the SND and DM2 data in the energy regiony^s 

G\ ■ up to 2 GeV. The u-meson parameters: = 782.79±0.08±0.09 MeV, 

F u = 8.68 ± 0.04 ± 0.15 MeV and a(uj -> 3tt) = 1615 ± 9 ± 57 nb were 
obtained. It was found that the experimental data cannot be described 
' by a sum of only uj, <f>, uj' and uj" resonances contributions. This can be 

\ interpreted as a manifestation of p — > 37r decay, suppressed by G-parity, 

with relative probability B[p -> 3tt) = (1.01 ±g;|| ±0.034) x 10~ 4 

">< 
Ph 



I. INTRODUCTION 



- 7r + 7r _ 7r° process in the energy region-^/s < 2 GeV is 
determined by the transitions of light vector mesons V (V = oj,4>,uj' ,uj") into the 
final state: V — > 7r + 7r~7r°. The V — > 7r + 7r~7r° branching ratios for vector mesons 
with isospin J = are large: B(ui -> 3tt) ~ 0.9, B(0 -> 3tt) ~ 0.15 1 , 5 (a/ -> 3tt) ~ 
1, — > 37r) ~ 0.5 2 , and thus the e + e~ — > 7r + 7r~7r° cross section measurements 
are important for study of these resonances. The pir intermediate state dominates 
in these transitions. The V — > 7r + 7r~7r° transition can also proceed via mechanisms 
suppressed by the G-parity: V — > unr° — > 7T + 7r~7r° 3,2 or V — > ,07r — > 7T + 7r~7r° 
(1/ = p,p',p"). The studies of the e + e~ — > 7r + 7r~7r° reaction allow to determine 
the vector mesons parameters and provide information on the OZI rule violation 
in <j) — > 37r decay and on the G-parity violation in the processes p, p'^"^ — ► 37r. 

The process e + e~ — * tt + tt~tt q in the energy region y/s below 2200 MeV was 
studied in several experiments during the last 30 years. The w-meson production 
region was studied in the Ref. 4-10 and studies of the <^-meson energy domain were 
reported in Ref. n ~ 16,8 . In Ref. 17-19 ' 8 the e + e~ — > ir + ir~ir° cross section was studied 
in the wide energy region-y/s from 660 up to 1100 MeV and in the Ref. 20 the upper 
limit was imposed on the G-parity suppressed decay p — > 3it. The e + e~ — > ir + ir~ir° 
cross section measurements in the uj' and uj" resonances energy region (/s = 1100 
- 2200 MeV) were reported in Ref. 21 " 24 ' 19 ' 8 ' 25 . 
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Recently the process e+e~ — > 7r + 7r~7r was also studied with the Spherical Neutral 
Detector (SND) 26 ^ 28,2 , the process dynamics was analyzed and the cross section 
was measured in the energy region y/s from 980 to 1380 MeV. Here we present the 
e+e~ — ► 7T + 7T~7r cross section measurement in the energy region-^/s below 980 MeV. 
The obtained cross section was analyzed together with the SND 27,2 and DM2 25 data 
in the energy region up to 2000 MeV. 

II. EXPERIMENT 

The SND detector 29 operated from 1995 to 2000 at the VEPP-2M 30 collider in 
the energy range y/s from 360 to 1400 MeV. The detector contains several subsys- 
tems. The tracking system includes two cylindrical drift chambers. The three-layer 
spherical electromagnetic calorimeter is based on Nal(Tl) crystals. The muon/veto 
system consists of plastic scintillation counters and two layers of streamer tubes. 
The calorimeter energy and angular resolutions depend on the photon energy as 
(te/E(%) = 4.2%/{/£(GeV) and a^ 9 = 0.82%/£(GeV) 0.63° The tracking 
system angular resolution is about 0.5° and 2° for azimuthal and polar angles re- 
spectively. 

In 1998 - 2000 the SND detector collected data in the energy regiom/i < 980 MeV 
with integrated luminosity about 10.0 pb _1 . For the luminosity measurements, the 
processes e + e~ — ► e+e~ and e+e~ — ► 77 were used. In this work the luminosity 
measured by using e + e~ — ► 77 was used for normalization, because in the p-meson 
energy region the contribution of the e + e~ — > tt + it~ background to the e + e~ — > 
e+e~ process is rather large. The systematic error of the integrated luminosity 
determination is estimated to be 2%. Since luminosity measurements by e + e~ — > 
e + e~ and e + e~ — > 77 reveal a systematic spread of about 1%, this was added 
to the statistical error of the luminosity determination in each energy point. The 
statistical accuracy was better than 1%. 

The beam energy was calculated from the magnetic field value in the bending 
magnets and revolution frequency of the collider. The relative accuracy of the 
energy setting for each energy point is about 0.1 MeV, while the common shift 
of the energy scale for all points within the scan can amount to 0.5 MeV. In the 
three energy points in the vicinity of the w-resonance peak the beam energy was 
calibrated using resonant depolarization method 31 . The accuracy of center of mass 
energy calibration is 0.04 MeV. In order to correct the calculated beam energy, the 
common shifts of the energy scale in the experimental scans were the free parameters 
in the analysis and varied relative to the calibrated energy points. The beam energy 
spread varies in the range from 0.08 MeV at^s = 440 MeV to 0.35 MeV at^s = 970 
MeV. 



III. DATA ANALYSIS 

A. Selection of e + e~ — > 7r + 7r~7r° events 

The data analysis and selection criteria used in this work are similar to those 
described in Rcf. 27,28,2 . During the experimental runs, the first-level trigger 29 selects 
events with energy deposition in the calorimeter more than 180 MeV and with two 
or more charged particles. During processing of the experimental data the event 
reconstruction is performed 29,27 . For further analysis, events containing two or more 
photons and two charged particles with \z\ < 10 cm and r < 1 cm were selected. 
Here z is the coordinate of the charged particle production point along the beam axis 
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(the longitudinal size of the interaction region depends on beam energy and varies 
from 2 to 2.5 cm); r is the distance between the charged particle track and the beam 
axis in the ? — <f> plane. Extra photons in e + e~ — > 7r + 7r~7r events can appear because 
of the overlap with the beam background or nuclear interactions of the charged 
pions in the calorimeter. Under these selection conditions, the background sources 
are e + e~ — > a>7r° — > 7r + 7r~7r 7r , e + e~ — > e + e~7, e + e~77, ir + ir~(j), /i + /i~(7) 
processes, cosmic and beam backgrounds. 

The polar angles of the charged particles were bounded by the criterion: 15° < 
6 < 165°. To suppress the cosmic and beam backgrounds, the following cuts were 
applied: E neu > 100 MeV, 2 < N y < 3, |Az| < 3 cm and ip > 20°. Here E neu is the 
energy deposition of the neutral particles, N y is the number of detected photons, 
Az = z\ — Z2, and z\, z 2 are z-coordinates of the charged particles tracks, tp is the 
angle between two charged particles tracks. 

To suppress the e + e~ — > e + e~77 events, an energy deposition of the charged 
particles in the calorimeter E c h a wa s required to be small: E c h a < 0.5 -\/s. 

To reject the background from the e + e~ — > 7r + 7r~(7), and e + e~7, the 

following cut was imposed: \A<j>\ > 5°. Here A<j) is an acollinearity angle of the 
charged particles in the azimuthal plane. 

For events left after these cuts, a kinematic fit was performed under the following 
constraints: the charged particles are assumed to be pions, the system has zero total 
momentum, the total energy is s/s, and the photons originate from the ir° — ► 77 
decays (Fig.l). The value of the likelihood function \\-k (Fig. 2) is calculated during 
the fit. In events with more than two photons, extra photons are considered as 
spurious ones and rejected. To do this, all possible subsets of two photons were 
inspected and the one, corresponding to the maximum likelihood was selected. The 
kinematic fits were also performed under assumptions that the e + e~ — > 7r + 7r~7, 
e + e~ — > /i + /i _ 7 or e + e~ — > e + e~7 events with extra photons were detected and the 
values of the corresponding likelihood functions xiv^i xi^7 an d X2e-y were calculated. 
After the kinematic fits, the following cuts were applied: xiL < 20, xl-w-y > 20, 
X2fi- ( > 20 and x\e~i > 20, the polar angle 6 1 of at least one of the photons, selected 
by the reconstruction program as originated from the tt° decay, should satisfy the 
following criterion: 36° < Q 1 < 144°. In the energy region yfs < 730 MeV, for 
additional suppression of the background, the cut p/E > 0.5 - 0.3 (aty/s — 720 - 
440 MeV) was applied. Here p and E are the charged pions momentum and energy 
calculated after the kinematic reconstruction. For additional suppression of the 
e+e~ — > U7T — ► 7r + 7T _ 7r 7r° background, the criterion N 7 — 2 was applied for the 
energies Vs > 900 MeV. 

The angular distributions of particles for the selected events are shown in 
Fig. 3, 4, 5, 6 and 7, while Fig. 8, 9, 10 and 11 demonstrate the photon and pion energy 
distributions for the same events. The experimental and simulated distributions are 
in agreement. 

B. Background subtraction 

The number of background events was estimated from the following formula: 

V fcfc3 ( S )-^a ffi ( S )e 2 ( S )/L( S ), (1) 

i 

where i is a process number, <Jm(s) is the cross section of the background process 
taking into account the radiative corrections, IL(s) is the integrated luminosity, 
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ej(s) is the detection probability for the background process obtained from simu- 
lation under selection criteria described above. For the e + e~ — > lo-k background 
estimation, the cross section obtained in SND experiments 32 was used. 

To estimate the accuracy of background events number determination, the Xz* 
distribution in the energy regiony / s > 870 MeV (Fig. 12) was studied. The experi- 
mental X3tt distribution in the range < < 40 was fitted by a sum of background 
and signal. The distribution for background events was taken from the simulation, 
while for the signal e + e~ — ► 7r + 7r~7r events - by using experimental data collected 
in the vicinity of the u> meson peak. As a result, the ratio between the number of 
background events obtained from the fit and the number calculated according to 
Eq.(l) was found to be 2.0 ± 1.2. The error was estimated by varying the selection 
criteria. Taking into account this ratio, the number of background events obtained 
from Eq.(l) was multiplied by a factor of 2 in all energy points and the accuracy of 
the determination of the number of background events was estimated to be about 
60%. 

The numbers of e + e~ — > n + Tr~ir° events (after the background subtraction) and 
background event numbers are shown in Table I. 

To check the accuracy of background subtraction in the energy region ^/s < 730 
MeV, the data were analyzed in a different way. The kinematic reconstruction 
was performed under the following constrains: the charged particles are pions, the 
system has zero total momentum, the total energy is The constrain that the 
photon originated from the 7r° — ► 77 decay was not used. In events with more 
than two photons, extra photons are considered as spurious ones and rejected. The 
value of the likelihood function X2ir2j 1S calculated during the fit. After the fit, the 
following cut was applied: < 20. For selected events the two-photon m 77 

invariant mass spectra (Fig. 13) were fitted by a sum of background and signal (m 77 
is calculated after the kinematic fit). The shape of the distribution for e + e~ — > 
7T + 7r~7r° events was obtained by using experimental data collected in the vicinity 
of the w-resonance. For the background the uniform distribution was used (other 
assumptions about the shape of the background spectrum do not change the fit 
results) . The cross sections obtained by using two different methods of background 
subtraction are in agreement (Tablell). The analogous check was performed for 
the energy region y/s > 900 MeV. The results of different approaches are again in 
agreement. 



C. Detection efficiency 

The detection efficiency of the e + e~ — > 7r + 7r~7r°(7) process was obtained from 
simulation. To take into account the overlap of the beam background with the signal 
events, background events (experimental events collected when the detector was 
triggered with an external generator) were mixed with the simulated events. The 
detection efficiency for events without 7-quantum radiation by the initial particles 
is about 0.35 (Table I). The detection efficiency dependence on the radiated photon 
energy is shown in Fig. 14. The efficiency decrease with the rise of the radiated 
photon energy is due to the selection criterion xjL < 20, which involves the energy 
and momentum conservations in the e+e - — > 7r + 7r _ 7r° process. 

Inaccuracies in the simulation of the distributions over some selection parameters 
lead to an error in the determination of average detection efficiency. To take into 
account these uncertainties, the detection efficiency was multiplied by correction 
coefficients, which were obtained in the following way 27 ' 2 . The experimental events 
were selected without any conditions on the parameter under study, using the se- 
lection parameters uncorrelated with the studied one. The same selection criteria 
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were applied to simulated events. Then the cut was applied to the parameter and 
the correction coefficient was calculated: 



n/N 

6 =^?M> (2) 

where N and M are the number of events in experiment and simulation respectively 
selected without any cuts on the parameter under study; n and to are the number of 
events in experiment and simulation when the cut on the parameter was applied. As 
a rule, the error in the coefficient S determination is connected with the uncertainty 
of background subtraction. This systematic error was estimated by varying other 
selection criteria. 

The inaccuracy in xi-n distribution simulation (Fig. 2) is the main source of 
uncertainty in the detection efficiency determination. The correction coefficient 
S x 2 = 0.95 ± 0.02 (Fig. 15) was obtained by using data collected in the vicinity of 
the u> resonance. The error due to uncertainty in simulation of other parameters is 
estimated to be 1.7%. In the energy region y/s > 900 MeV the additional selection 
criterion 7V 7 = 2 was applied. The uncertainty due to this cut was estimated to 
be 3%. The systematic error of the detection efficiency determination is 2.7% in 
the energy region y/s < 900 MeV and is about 4.1% at-^/i > 900. The detection 
efficiency after the applied corrections is shown in Table I. 



IV. THEORETICAL FRAMEWORK 

In the framework of the vector meson dominance model, the cross section of the 
e+e~ — ► 7r + 7r~7T° process is 

da 4na \p + x p_| 2 2 

-m m+ • \F\ Z , (3) 



dniodm+ s 3 / 2 127r 2 i/s 

where p+ and p- are the 7r + and 7r~ momenta, too and to + are 7r + 7r~ and 7r + 7r° 
pairs invariant masses. The form factor F of the 7* — > 7r + 7r~7r° transition has the 
form 



IFI 2 = 



A 



2 

(4) 



The first term in Eq.(4) takes into account the 7* — > pn — > tv + tt tt° transition 
(Fig. 16 a), which dominates in the process under study 28 . Here 



D p {mi) = to 2 , - to- - imiT pi (mi), T p i(rrn) = I I T pi 
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m i I 9 \ Qi ( m p i ) 



<?o(™ 2 ) = i(TO 2 -4m 2 ) 1 / 2 , q±(m 2 ) = ^- [(TO 2 -(m 7r o+TO 7r ) 2 )(TO 2 -(TO 7r o-TO 7r ) 2 )] 1/2 , 



m. 



= \Js + to 2 + 2to 2 - ml- m\, 



where to_ is the n tt° pair invariant mass, m n o and m n are the neutral and charged 
pion masses, i denotes the sign of the p-meson {-kit pair) charge. The p a — ► -k + it~ 
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and p^ — > ti^tt transition coupling constants could be determined in the following 
way: 



2 



67rTO 2 o r p o 



6nm 2 p± r p ± 
q±(m p ±) 3 



Experimental data 28 do not contradict to the equality of the coupling constants 
q 2 n0 „„ = Q 2 ± ■ In this case the p° and p^ meson widths are related as follows: 



m p0 q±(m p ±y 



m 2 ± q {m p of 



(5) 



In the subsequent analysis we assume that g p o n7r 



taken from the SND measurements 28 : T p o 



o 2 + , and the width values were 



149.8 MeV, T p ± 



150.9 MeV. The 

neutral and charged p mesons masses were assumed to be equal and were also taken 
from the SND measurements 28 m p = 775.0 MeV. A factor Z(m) = 1 — isi$(m, s) 
takes into account the interaction of the p and ir mesons in the final state 33 (Fig. 16 
d), parameter Si = 1±0.2 corresponds to the prediction of Ref. 33 , where the concrete 
form of the <I>(to, s) function can be found. In experimental studies of the tttt mass 
spectra in the e + e~ — > 3tt process at^/s ~ m^ 28 we had obtained si = 0.3±0.3±0.3. 
This result is consistent with zero, but also does not contradict to the prediction of 
Ref. 33 . 

The second term in Eq.(4) takes into account the possible transition 7* — > 
p'(")n — > 7r + 7T~7r° (Fig. 16 c). This term can be written as A p7r (s) ■ a 37r , where 



037T 



V">-( s ) 



.9 



y . 

^ r> . ( „)(m 4 ) 



(s) ^ D 



In the analysis, the a 3Tr amplitude was assumed to be a real constant. From the rnr 
mass spectra analysis, in the process e + e~ — ► 37r at^/s ~ m^ 28 , it was found that 
a 37r = (0.01 ± 0.23 ± 0.25) x 10~ 5 MeV" 2 . 

The third term in Eq.(4) takes into account the 7* — > unr° — > ir + ir~ir a transition 3 
(Fig. 16 b). The polarization operator of the p — u) mixing satisfies lm(U p(Aj ) <C 
Re(n pw ) 33 - 34 , where 



Re(n pw ) 



T p o (m w ) 



B(ui — > 7T+7T ) 



(m 2 - to 2 ) - «TO w (r w - r p o(m w )) 



, (6) 



so we have assumed Im(n pu) ) = in the subsequent analysis. 

The e+e~ — > 7r + 7r _ 7r° process cross section can be written in the following form: 



(7) 



where 



4?ra w c ^ 



A P tt(s) 



(8) 



- 47ra w r \ 

&u>7r— >3tt — ^3/2 ^^v^V 



(9) 



f int — 



372 



^(5)^(5)^(5) + A*(s)A un (s)W* nt (s) 



(10) 
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The phase space factors W p7r (s), W W7r (s) and Wj„ t (s) were calculated as follows: 

s/s-m^a m+ ax ('"o) 



127r 2 yi 



J m dm a J m + dm + \p + x p_\ 2 ■ ^ 

2m. m ~»(m ) *=+.0: 



i=ffl t - D p (mi)Z(mi) 



(11) 



v/s-m >0 m™ ax (m ) 

m dm n 



2nit m~"(m ) 



D p (m a )D u (mo) 



(12) 



Wi„ t («) 



/s-m i0 m™ ax (m ) 

///or/'/'. ' '" 7 ' ' ' ' 



12^/5 / m ° rfm ° / 

2^tt m r T' zn (rr< 



m+dm+\p+ x p_ 



n 



m~"(mo) 



£)p(mo)£> w (m ) 



.9p* 7T7T 



i=+,0, 



■ p (mi)Z(mi) 



+ a 37l 



(13) 



Amplitudes of the 7* — > /97r transition have the form 

1 r F m^m y< 7(V -> 3tt) e^ y C^(i,r ) 



^4p-7r(s) 



7ra 



E 



V— LU,p,<p,U}' ,Uj'' 



D v (s) 



(14) 



where 



£V(s) - ™ 2 y - s - i>/nv(*), r v («) = ^r(v -> /, *). 

/ 

Here / denotes the final state of the vector meson V decay, my is the vector meson 
mass, IV — IV(my). The <j> meson mass and width were taken from the SND 
measurements = 1019.42 MeV, T = 4.21 MeV 27 . The following forms of the 
energy dependence of the vector mesons total widths were used 



r P («) 



™ P qo(s) 3 
s q {m p ) 



12tt 



1 1 ' = —7^4^5(0; 7r+7r-)C^(5,r )+ ^^' r.B^ - 7r° 7 )C7 2 ( S , r )+ 



s qvim^f 



r*(«) 



s q K ±(m 4 ,y 



-T^^ K+K-)Cl K+K -(s,r )+ 



m l q K o{s) 3 
s q K a( m4 ,y 



rT^B{<t> -+ K s K L )Cl KsKL (s, r ) + 
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2 ln „ W pA s ) 



r w /(s) = i\yC w , (s,r- ) 



r W "(s) = r^c2 ( s ,r ) 



Here is a coupling constant of p — > W7r° transition, q W7r , <?k±, g^o, q 7r7 and 
the w meson, kaon, 77 meson and pion momenta, W wir7r (s) is the phase 
space factor of the w7T7t final state 35 , Cypp(s, ro) and Cvvp(s,tq) are the form 
factors which restrict too fast growth with energy of the partial widths, so that 
y / ir(s) — > cons/j as s — > 00. According to Ref. 36 these form factors can be written 
as follows 



Cvy{p)p( s , r o) — 



1 + (r m v ) 2 
1 + (»V5) 2 ' 



Cvp P (s,r ) =4 



1- + (^ogp(^y)) 2 
l + (r qp(s)) 2 ' 



(15) 



where gp is the momentum of the pseudoscalar meson, ro is the range parameter 
(its value was taken to be the same for all decays). 

The relative probabilities of the decays were calculated as follows 



a(V) 



x 



ml 



In particular: 



B(w -» X) 



a(uj -» X) 
ct(w) 



ct(w) = 



a(u> — > 37r) + cr(ti) — > 7r°7) 

1 - S(W 7T+7T-) ' 



In further analysis we have used a(ui — > 7r°7) = 155.8 nb, cr(<^) — » i^+if - ) = 1968 
nb, ct(0 — > K S K L ) = 1451 nb and er(</> — » 777) = 54.8 nb obtained in the SND 
experiments 37 ' 27,38 . 

(p^v is a relative interference phase between the vector meson V and w, so (p^ = 
0°. Phases 4>^v can deviate from —180° or 180° and their values can be energy 
dependent due to mixing between vector mesons. For example, the phase <j) U( f, 
was found to be close to 180° 27 ' 2 in agreement with the prediction 39 <fiu<f> = $(s) 
($(m^) ~ 163°), where the function $(s) is defined in Ref. 39 . In Ref. 2 it was shown 
that (\> ULO i ~ 180° and </> wa) » ^ 0°, so in this work these two phases were fixed on 
those values. 

Taking into account the p — u) mixing, the w — > pir and p ^ pn transition 
amplitudes can be written in the following way 3,34 



where 

e(s) = 



-n„ 



DUs)-Dp(s) 



(0) o (0) 



\gvj\ 



(0) 



+ 



-(0) (0) r (0) 



(0) 



s(s) 



3m 3 v T v B(V-> e+e") 



1 1/2 



4ira 



\gv P 7r\ 



(16) 



47Tl>g(y -> /77T) 



1/2 



The superscript (0) denotes the coupling constants of the pure, unmixed state. The 
Eq.(16) can be rewritten as follows 



A-i^j — > pn 



A 



1 



E 



Yym\^m v o(V -> 3tt) fy p ^{s)C v p7r {s, r ) 



I>v(«) 
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where 



fvpn(s) 



rv P ir{m v )' 



(0) 

^(S) = l + ^£(s)~l + 



m p i T(p —> e+e ) 
m^T(uj — > e+e~) 



1/2 



e(s), 



(o) 

r W7r (s) = ^--e{s) ~ -e(s) 

.9o)p7T 



If the p — > 37r transition proceeds only via p — u mixing, that is g pp \ = 0, then 
<j) up ~ —90° and almost does not depend on energy, besides 



a{p — > 37r) ~ cr((jj — » 37r) 



m 2 W plx {m p ) r pplT (m p ) B(p — > e+e ) 
T p W^m^) r W(97r (m w ) B(w -> e+e") ' 



For the 7* — » u>7r° transition amplitude, the model which gives satisfactory de- 
scription of the relative phase between it and the A pn (s), Eq.(14), amplitude 2 was 
used: 



m 3 p T p B(p-^ e+e )g puJ7I r y m 2 v c^Q) e^ v 

D p(s) + v£/,p" Dv{s) VtiMw) 



where <j) pp = 0°, m p « = 1480 MeV, T p , = 790 MeV, a(p' -» wtt°) = 86 nb, </y = 
180°, rry- = 1640 MeV, I> = 1290 MeV, -> wtt°) = 48 nb, = 0° 

5pw7r = 16,8 GeV -1 and 



r P '(») (s) = r 



V. CROSS SECTION MEASUREMENT 

From the data in Table I, the cross section of the process e+e~ — > 7r+7r~7r° can 
be calculated as follows: 

- ilwfw • (18) 

where N 3n (s) is the number of selected e+e - — > 7r+7r~7r°(7) events, IL(s) is the 
integrated luminosity, £(s) is the function which takes into account the detection 
efficiency and radiative corrections for the initial state radiation: 

E raax 

J cr 37r (s, E 1 )F{s, E y )e(s, £" 7 )dE 7 

m = -5 —7-, • (19) 

Here £ 7 is the emitted photon energy, F(s, E 7 ) is the electron "radiator" function 40 , 
e(s, E^) is the detection efficiency of the process e+e - — > 7r+7r _ 7r°(7 rad ) as a func- 
tion of the emitted photon energy and the total energy in the e+e~ center of mass 



9 



system, <r 37T (s) is the theoretical energy dependence of the cross section given by 
the equation (7). 

To obtain the values of £(s) at each energy point, the visible cross section of the 
process e + e~ — ► Tr + -K~Tr a (j ra d) 

was fitted by theoretical energy dependence 

The following logarithmic likelihood function was minimized: 

(erf s - of 1 ) 2 



* 2 = E- 



A? 



where i is the energy point number, Aj is the error of the visible cross section a ms . 

In the fit the </>, w', u>" meson parameters (the mass, width, branching ratios 
of main decays ) was fixed at their values obtained in the SND experiments 27 ' 2 , 
other parameters were fixed as follows ro = 0, a 37T = 0, Si = and (j)^^ = $(s) 
($(m^) ~ 163°). The Eq.(14) was written in the following form: 



where is a complex constant. 

The u{uj — > 37r), r w and m u were the free parameters of the fit. The £(s) values 
were obtained from the approximation of the experimental data in several models: 

1. a{p-> 3tt) = 0, C 37r = 

2. er(/9 — » 37r) and W(9 are free parameters, = 

3. C-i-K is a free parameter, a(p — > 37r) = 

4. <t(p — ► 3tt), and C3,r are free parameters 

The fits were also performed under the same assumptions, but with o(uj' — > 37r) = 
and er(w" -> 3tt) = 0. 

The values of £(s) actually do not depend on the applied model. The largest 
model dependence, about 1.5 — 2%, was found at-^/s from 800 to 840 MeV. Using 
the obtained £(s) values, the cross section of the e + e~ — ► 7r + 7r~7r process was 
calculated (Table III ). The systematic error of the cross section determination at 
each energy point -^/i is equal to 

<J S ys — CTeff © &IL © Vmod{s) <Jbkg(s)- 

Here a eff = 2.7% at^ < 900 MeV and <r eff = 4.1% at-s/s > 900 MeV, a IL = 2%. 
They are systematic uncertainties in the detection efficiency and integrated lumi- 
nosity, which are common for all energy points. The model uncertainty <J mo d(s) 
was obtained from the difference of £(s) values determined for the models men- 
tioned above. The error <Jbk g {s) takes into account the inaccuracy (~ 60%) of the 
background subtraction and depends on the beam energy. 
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VI. THE e+e" -f tt+tt-tt CROSS SECTION ANALYSIS. 



The cross section measured in this work (Table III) was analyzed together with 
the e + e~ — > 7r + 7r~7r cross section measured by SND in the energy region-^/s from 
980 up to 1380 MeV 27 ' 2 and with the DM2 results of the e+e~ -» tt+tt^tt and 
e + e~ — » W7r + 7r _ cross sections measurements in the energy region-^i from 1340 to 
2200 MeV 25 . 

The e + e~ — > 7r + 7r~7r° cross section was fitted by the expression (7). The e + e~ — > 
o>7r + 7r~ process cross section was written in the following way: 



1 



r w »m 2 „ vM^" W7r+7r-)m w » /W a;7r7r (s)Cyyp(s,ro) 



(21) 



Here we have neglected the ui' — > o>7T7r contribution. To calculate B(lo" — > Loirir), 
we have assumed a(u" — > o>7T7r) = 1.5 • o-(u/' — > W7r + 7r~). 

The cross sections of the e+e - — > 7r + 7r~7r° and cj7r + 7r~ processes, measured by 
SND and DM2, were fitted together. The function to be minimized was 



where 



2 2 2 2 

Xtot — Xz-n(SND) + X3tt(_DM2) + ^ > nr(OM2) ' 



2 _ V- f a il ND) ( s ) - VSit (s) 

X3tt(SND) — 2^ I 



4° ND) (s) 



2 _ vY C-g^ M2) (g) ~cr 37r (s) 

X^(DM2) _ Z^I .(DM2), x 

s \ A 3tt l s J 



/ (DM2), x / \\ 

X W 7r7r(£)M2) — I .(DM2), s I 

Here v\ S ^^Tr) M2 ^ ( s ) are ^ ne experimental cross sections, A are their uncertainties, 
C is a coefficient which take into account the relative systematic bias between 
SND and DM2 data. The errors A 37r ' include both the statistical <J s tat and the 

systematic errors: /S^ ND ^ = a s tat © &mod © o'bkg- I n Ref. 2 the C coefficient was 
estimated to be C = 1.54. In the analysis that follows we have fixed this coefficient 
at 1 or 1.54. 

In the fittings m w , T w , a(u> — > 37r), <t(<^) — > 3-7r), </> W 0, m^/, IV, cr(u/ — > 37r), 
Tow", r^", cr(u/' — » 37r) and er(a/' — > o>7r + 7r~) were the free parameters. The 
approximations were performed under the following assumptions about the phase 
space factor for the 7T + 7r _ 7r° final state: 

1. s 1 = 0, a 37T = 0; 

2. si = 1, a 3Tr = 0; 

3. si = 0, a 37r = -4 x 10~ 6 MeV" 2 ; 

4. si =0,a 3 „=4x 10~ 6 MeV" 2 . 
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The nonzero values of the a 3T amplitude are the upper limits imposed on the 90% 
confidence level by using the SND result reported in Ref. 28 . The approximations 
were also performed without taking into account the contribution from the e + e~ — > 
con — > 7T + 7r _ 7r° process, i.e. by assuming = a p7T . The difference in the fit 
results was included in the model uncertainty. 

The fittings were performed under the following model parameters: 

1. r = 0, <r(p -> 3tt) = 0; 

2. ro = 0, <j) up and cr(p — > 37r) are free parameters; 

3. a(p —> 3it) = 0, ro is a free parameter. 

The results of the fits are shown in Table IV and in Fig.17, 18,19,20,21,22. In 
the Table IV, xl, x\ S nd) (xL(swd) = xl + x\ S nd)) and denote the x 2 
values for the energy rcgions-^/s below and above 970 MeV and in the energy range 
880 <y/s < 970 MeV. The x 2 values in the first model (column 1 in the Table IV) is 
too large and this model contradicts to the experimental data. The second and third 
models (columns 2 and 3 in Table IV) are in agreement with the experimental data. 
The Xtot vahie for the second model is less than for the third one. If <r(p — ► 3tt), 
(j)(jjp and r are the free parameters in a fit, then they are found to be equal to 
a{p -» 3tt) = 0.11±g;g| nb, wp = -136±}g degree and r = 0.2 ± 0.3 GeV" 1 . 
In this case the value of the range parameter ro turns out to be rather small and 
consistent with zero. The value of Xsnd f° r the second model is less than for the 
third one, and vice versa the xl f° r the third model is less than for the second. In 
the energy range 880 <-y/i < 970 MeV (Fig. 20) the fitted curves for the second and 
third models exceed the experimental points in average by about 1.5 and 1 errors 
respectively. The difference between models 2 and 3 is also seen in the energy 
regions^ < 720 MeV and^i > 1100 MeV (Fig.17 and 22). The xl wn m M 2) for the 
third model increases by a factor of 2 in comparison with the second model. The 
results of the e+e~ — ► W7r + 7r~ cross section fits are shown in Fig. 23. In case when 
the form factors (15) are used, the theoretical curve poorly describes experimental 
points at the left slope of the resonance . The CMD2 results of the e + e~ — ► ujtt + tt^ 
reaction studies 41 are also presented in Fig. 23. These data agree better with the 
second model. 

If the relative bias between SND and DM2 measurements is not assumed, the 
*L(DM2) valuc is rather lar £ e: xL(DAf2)/ Ar /« = ( 37 ^ 40)/18. Here N flt is 
the number of fitted experimental points (Table IV). A rather large scale fac- 
tor C — 1.54 is required to concert the SND and DM2 data, and in this case 
X 2 . n (DM2)/Nfit = (22-=-27)/18. In order not to guess about relative systematics be- 
tween the SND and DM2 experiments, the fits described above were redone assum- 
ing C = 1, but without taking into account e + e~ — > 7r + 7r~7r° cross section measured 
in DM2 experiments. (Table V). The parameters m w », r<y and a(uu" — > aj7r + 7r~) 
were obtained from the fitting to the cross section e+e~ — > ljtt + it~ reported by 
DM2, and m w >, T u i, <j(uj' — > 3ir), o{uj" — > 37r) were obtained by using SND data 
only. In this case the first model (column 1, Table V) agrees with experimental 
data, but agreement is significantly better if the fits are performed in the second 
or third models (columns 2 and 3 in Table V). The Xtot value for the third model 
is slightly bigger than for the second one. In this approach the fitted curve is in 
conflict with the DM2 measurenment of the e+e~ — > 7r + 7r~7r cross section (Fig. 24). 
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VII. DISCUSSION 



Comparison of the e+e — > 7r+7r 7r° cross section obtained in SND experiments 
with other rcsu lts 17 ' 8, 19,10,15,16 is shown in Fig.25,26,27 and 28. The DM1 results 17 
are in agreement with the SND measurements. The ND results 8 ' 19 agree with SND 
data in the energy region yfs < 930 MeV, while fory^ > 930 MeV ND points lay 
on about two errors lower than SND ones. In the vicinity of the u> resonance peak 
()/s ~ 780 MeV) the SND cross section exceeds the CMD2 measurements 10 , while 
in the <j> meson energy region the SND and CMD2 results 15,16 are in agreement. 

The u) meson parameters m u , T w , a(uj — » 37r) were measured through study of 
the e + e~ — > Tr + ir~ir° cross section. The ui meson mass was found to be 

m u = 782.79 ±0.08 ±0.09 MeV. 

Here the systematic error is related to the accuracy of the VEPP-2M energy scale 
calibration by resonant depolarization method, 0.04 MeV, and to the model un- 
certainty, 0.08 MeV. The SND measurement in comparison with the results of 
experiments 10 ' 42,43,6 ' 44,17 and the world average value m^ 1 is shown in Fig. 29. SND 
result is in agreement with the CMD2 measurement 10 , and differs from world av- 
erage by about 1.3 standard deviations. The maximum difference, about 3.4 stan- 
dard deviations is between the SND result and the Crystal Barrel measurement 
= 781.96 ±0.21 43 . 
The following value of the w meson width was obtained: 

r w = 8.68 ±0.04 ±0.15 MeV. 

The systematic error is related to the model dependence and to the accuracy of 
energy determination. The comparison of this value with the results obtained in 
Ref. 10 ' 42 ' 7 ' 9 ' 6,17,5 and with PDG world average value 1 is shown in Fig.30. The SND 
result agrees with other measurements. 

The parameter a(uj — > 37r) was found to be 

a(u -> 3tt) = 1615 ± 9 ± 57 nb. 

The systematic error includes the systematic uncertainties in the detection efficiency 
and luminosity determinations, 55 nb in total, and the model dependence 13 nb. 
The comparison of the obtained value with other experimental results 10,7,9 ' 6 ' 17 ' 5 and 
with PDG world average 1 is shown in Fig. 31. The SND result exceeds the central 
values of the majority of the previous measurements. It differs by less than one 
standard deviation from the results in Ref. 7 ' 9 ' 5 , by about 1.4 standard deviations 
from the DM1 measurement 17 , by 2 standard deviations from the OLYA result 6 
and PDG world average <t(ui — > 37r) = 1484 ± 29 nb. The difference from the most 
precise measurements, done by CMD2 10 and SND, is about 2.5 standard deviations. 

Using the SND result a{uo -> tt°7) = 155.8±2.7±4.8 nb 37 , the ratio of the partial 
widths of the uo — > 7r°7 and uj — > 3tt decays was calculated 

T((jJ -> 7T°7) 

— ^ = 0-097 ± 0.002 ± 0.005 

T{uj — > 37r) 

This value agrees with PDG world average 1 and with other experimental 
results 45 ' 7 ' 44 ' 5,46 (Fig.32). 

Using a(co — > 3w), measured in this work, the SND result of the u> — > 7r°7 decay 
study 37 and PDG world average value B(u> -» tt+tt") = 0.0170±0.0028 1 , the partial 
width of the u — » e+e~ decay and the u meson main decays branching ratios were 
obtained: 

T(lu -> e+e") = 0.653 ± 0.003 ± 0.021 keV, 
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B{uj -> e+e") = (7.52 ± 0.04 ± 0.24) x 10~ 5 , 
B(lo -» 3tt) = 0.8965 ± 0.0016 ± 0.0048, 

B(w -» tt°7) = 0.0865 ± 0.0016 ± 0.0042. 

Comparison of these results with PDG data 1 is presented in Table .VI. The value of 
B(u) — ► e+e~), calculated by using SND data, exceeds the world average by about 
2 standard deviations (by 8%). 

The e + e~ — ► 7r + 7r~7r° and e+e~ — » a)7r + 7r~ cross sections analysis show that 
the data cannot be described by a sum of u, <fi mesons and two u)', uj" resonances 
(model 1). The data can be satisfactory described with the model 3, which takes 
into account form factors (15), with constrained partial widths growth with energy. 
The range parameter of this form factors was found to be 

r = 2.5 ±J;| ±0.5 GeV -1 . 

This agrees with expected vector meson effective "radius": 2.5-3 GeV -1 47 . The 
second error is due to model dependence. The model 2, which takes into account 
the 7* — > p —f 37r transition, also satisfactory describes the experimental data. For 
parameters of this model, the following values were obtained: 

a{p -> 3tt) = 0.112 ±g$g ±0.038 nb, 

<p up = -135 ±H ±9 degree. 

Here the systematic error is due to model uncertainty. The above given <r(p — > 37r) 
value corresponds to the branching ratio B(p — > 37r) = (1.01 ±0 3! ±0.34) x 10~ 4 . 
Assuming that the p — > 37r transition proceeds via the p~uj mixing mechanism, the 
following values of the 7* — > p — > 37r process parameters are expected: ^ W(9 ~ —90° 
and cr(p — > 37r) = 0.05 - 0.07 nb. The a(p — » 37r) value obtained in the analysis 
agrees with the expected one, while differs from the expected value by about 
two standard deviations. 

In general, the model 2 seems to be more preferable than the model 3 due to the 
following considerations. The full data set for the e + e~ — ► 7r + 7r~7r° and ujn + ir~ 
cross sections is in somewhat better agreement with the second model. Approxi- 
mation of the e + e~ — > W7r + 7r~ cross section with the third model is poor (Fig. 23). 
The phase value, obtained by the fit in the second model agrees with the the- 
oretical prediction = $(m^) w 160° 39 , while the phase 4>oj<p, obtained by using 
the third model, exceeds the expected value (Fig. 34). But, unfortunately, available 
experimental data are insufficient to make a strict conclusion about observation of 
the p — > 37T decay. 

The parameter a{<p — > 37r) was found to be 

a((j) -» 3tt) = 657 ± 10 ± 37 nb. 

The systematic error includes the systematic uncertainties in the detection effi- 
ciency and luminosity determinations, 33 nb in total, and the model dependence 
17 nb. This value agrees with the results of our previous analysis 27,2 . The SND 
result agrees also with other measurements 16 ' 15 ' 8 ' 18 ' 17,14,12,11 and with PDG world 
average 48 (Fig.33). 

The fit within the models 2 and 1 (Tables IV and V) gave the result 

= 163 ± 3 ± 6 degree. 
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The systematic error is related to model dependence. The obtained result is in 
agreement with the theoretical prediction (Fig. 34) (j)^^ = W(s), ^(m^) = 163°, 
which takes into account the (j> — lo mixing 39 . The fits with the model 3 (column 3 
in Tables IV and V) gave the result cj)^^ = 190 ± 5 ± 10 degree, which exceeds the 
theoretical prediction. 

The conventional view on the OZI suppressed <j) ~ * 7r + 7r~7r° decay is that it 
proceeds through <f> — ui mixing, i.e. in the wave function of the 0-meson, which is 
dominated by s quarks, there is an admixture of u and d quarks: 

|0 >w |^ (0) > +e 4>UJ \^ 0) >, |0 (O) >= 88, M° } >= (uu + dd)/V2, 

£<t>u> ~ 0.05 is (f>— lo mixing parameter. An alternative to the <f>— lo mixing is the direct 
decay. In Ref. 49 ' 36 ' 34 ' 50 ' 51 it was shown that there are no serious reasons to prefer 
the 4> — uj mixing to the direct transition, and methods of determination of the <f> — > 
7r + 7!~~7r° decay mechanism were suggested. In particular it was proposed in Rcf. 51 
to analyze the T((j) — > e + e~) /T{lo — > e + e~) ratio. In this work the B(u> — » e + e~) 
based mainly on the SND data was obtained, and in the Ref. 27 the B(<f> — » e + e~) 
was measured by SND. We performed the analysis of the <fi and lo mesons lepton 
widths ratio based on the SND data only. To improve the accuracy of B(<f> — ► e+e~) 
determination, the SND results of the (f> — > decay studies 52 were used. The 

average of these measurements y/B(4> — > e + e~)B{<j> — > = (2.93±0.11) x 10~ 4 

agrees with B{4> -> e+e") = (2.93 ± 0.15) x 10~ 4 27 . Assuming B(<p -> e+e") = 
B(<p -» ^ + aO, onc g cts s O -» e+e") = (2.93 ± 0.09) x 10~ 4 . The ratio of the 
leptonic widths is equal to 

T(6 -> e+e") 
^e+ e - = ^ ^4 = 1-89 ± 0.08 

On the other hand this ratio can be written in the following form: 



Re+c 




(0) , (0) 



(0) _ (o) 
9~/u> £0o)<7 7 



(22) 



Using the obtained R e + e - value, equation (22) and the nonrelativistic quark model 
prediction: 

the <f> — lo mixing parameter w 0.06 was obtained. On the other hand, taking 
into account the equation g^ = + e^g^, and assuming g^ = and 

g^pn = 9u P tti we found w 0.06. Here g wpir and t/0 p7r coupling constants were 
calculated by using the SND results obtained in this work and in Ref. 27 , the phase 
space factors W p „.(m w ) and W pw (m^) were calculated assuming s± — and 03^ = 0. 
In this case the SND data agree with the 4>— lo mixing dominance in the <f> — * 7r + 7r~7r° 
decay. 

The ratio ffi/ffi = W% is valid if V>(0, m v ) wave function of the qq bound state 

at the origin behaves like | -0 (0, my) | 2 oc m v , that is corresponds to the Coulomb- 
like nonrelativistic potential. But experimental data on the vector mesons p, lo, 
(j), J/ip, T(15) leptonic widths support the \tp(0,my)\ 2 oc m v behavior. Indeed, 
according to Ref. 53 : 



167ra 2 



T(V -+ e+e-) = —^C v W0, m v )\< 

Illy 
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where Cy is the mean electric charge of the valence quarks inside the vector meson 
V (C 2 = 1/2, Cl = 1/18, C 2 = 1/9, C 2 ^ = 4/9, C 2 (1S) - 1/9). In the case 
of \ip(0,my)\ 2 oc m v , and in the absence of the mixing, the following ratios are 
expected: 

T(p -» e+e") : T(u> -» e+e") : r(<£ -» e+e") : r(J/V> -» e+e") : r(T(lS) -» e+e") = 

= C 2 p :(%:(%: C 2 ^ : C 2 (1S) = 4.5 : 0.5 : 1 : 4 : 1 (23) 

Using the SND r(w — > e+e~) and T(cp — » e+e~) values and the world average for 
the T(p -> e+e"), r(J/r/> -» e+e") and r(T(lS) -» e+e") 1 , we have found: 

r(p -» e+e") : r(w -» e+e") : T(> -» e+e") : r(J/V> -» e+e") : r(T(lS) -» e+e") = 

= 5.2 ± 0.2 : 0.495 ± 0.025 : 0.93 ± 0.05 : 3.98 ± 0.032 : 1, (24) 

This ratios agree with the expected (23). If |^»(0, my)\ 2 cx my, then /i>°V/^ = 
-^j2m (Aj /m ( f > and by using Eq.(22) w 0.015 can be obtained. In this case the 
coupling constant of the direct (f> — ► 7r + 7r~7r° decay 5^ w 0.7 • 50 p7r is required to 
describe the experimental value of B((p — > 37r), indicating the direct transition as 
the main mechanism of the decay. 

The following u' parameters were obtained from the fits (Table IV and V): 

m u > = 1400 ± 50 ± 130 MeV, 
r w < = 870 ±|gg ±450 MeV, 

<j{d -» 3tt) = 4.9 ± 1.0 ± 1.6 nb. 

The lo' decays mostly into 7r + 7r~7r° and its electronic width is T(u' — > e + e~) <~ 570 
cV. The w" parameters were found to be 

m ull = 1770 ±50 ±60 MeV, 
I> = 490 ±?gg ±130 MeV, 
-» 3tt) = 5.4± 2 ;° ±3.9 nb, 

o-(w" ^wtt+tt") = 1.9 ±0.4 ±0.6 nb. 

The to" resonance decays with approximately equal probability into 7r + 7r _ 7r° and 
lotttt: B(ui" — ► 37r) <~ 0.65, — > o)7T7r) ~ 0.35 and it has the electronic width 
T(ui" — > e+e - ) <~ 860 eV. The second errors shown above are due to the model 
uncertainty and possible bias between the SND and DM2 measurements. The lo' 
and lo" parameters, obtained in this work, are somewhat different from those ob- 
tained in our previous analysis 2 . In particular, the lo' and lo" full widths values 
decrease a little. This difference is attributed to the fact that a new data for the 
e + e~ — > 7r + 7r _ 7r° cross section below 1 GeV were added in the fits. Of cause, the 
obtained values are not precise measurements, they should be considered as rather 
approximate estimation of the lo' and lo" resonances main parameters. 

In the energy region 880 < \fs < 970 MeV, experimental points deviate from 
fitting curves (Fig. 20). The difference can be attributed to inadequacy of the applied 
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theoretical models, uncertainty of the u>' and w" resonances contributions. Maybe 
more accurate consideration of the vector mesons mixing is required. 
Analysis of the 

Using the e+e~ — > 7r + 7r _ 7r° cross section, obtained with SND detector in this work 
and in Ref. 27 ' 2 , the contribution to the anomalous magnetic moment of the muon, 
due to the ir + ir~ir° intermediate state in the vacuum polarization, was calculated 
via dispersion integral: 

«,.(Wi < 1-38GCV) = (^f) ^^ds, 
where s max — 1.38 GeV, s m i n = m v o + 2m n , K{s) is the QED kernel, 

Here cr(e + e~ — > 7r + 7r~7r°) is the experimental cross section, A/(s) and A/ l (s) are 
corrections due to leptonic and hadronic vacuum polarizations. The A/(s) was 
calculated according to Ref. 54 and A^(s) was obtained by using e+e~ — > hadrons 
total cross section. 

The integral was evaluated by using the trapezoidal rule. To take into account the 
numerical integration errors, the correction method suggested in Ref. 55 was applied. 
As a result we obtained: 

a M (37r,Vi < 1.38GeV) = (458 ± 2 ± 17) x 10 -11 . 

At present in BINP (Novosibirsk) the VEPP-2000 collider with the energy range 
from 0.36 to 2 GeV and luminosity up to 10 32 cm _2 s _1 (aty/s <~ 2 GeV) is under 
construction 56 . The e + e~ — > tt + -k~tt process studies in the energy region-^/s < 2 
GeV will be continued in future experiments with SND detector at this new facility. 



VIII. CONCLUSION 



The cross section of the process e + e~ — ► 7r + 7r~7r was measured in the SND 
experiment at the VEPP-2M collider in the energy region^ below 980 MeV. The 
measured cross section was analyzed in the framework of the generalized vector 
meson dominance model together with the e + e~ — > 7r + 7r~7r° and utir + ir~ cross 
sections obtained by SND and DM2 in the energy region 980 < \fs < 2000 MeV. 
The w-meson parameters: m w = 782.79 ± 0.08 ± 0.09 MeV, T w = 8.68 ± 0.04 ± 0.15 
MeV and <j(uj — » 37r) = 1615 ± 9 ± 57 nb were obtained. 

It was found that the experimental data cannot be described by a sum of ui, 
(f>, u)' and u)" resonances contributions. This can be interpreted as manifestation 
of p — > 37r decay suppressed by G-parity, with relative probability B(p — > 37r) = 
(1.01 ±[j;36 ±0.034) x 10 -4 . The relative interference phase between the w and p 
mesons was found to be equal to: 4>^ p = —135 ±\l ±9 degree. These parameters of 
the p — > 3tt decay are in agreement with the theoretical values expected from the 
p — lo mixing. 

Analysis of the T((f> — ► e + e~ )/r(oj — » e+e~) ratio and t/0 p7r and g wpT coupling con- 
stants obtained in SND experiments indicates that the direct transition is preferable 
to <p — lu mixing as the main mechanism of the 4> ~ * 7r + 7r~7r° decay. 

Using the e + e~ — ► 7r + 7r~7r° cross section, obtained with SND detector the contri- 
bution to the anomalous magnetic moment of the muon, due to the 7r + 7r _ 7r° inter- 
mediate state in the vacuum polarization, was calculated: a A1 (37r,- s /s < 1.38GeV) = 
(458±2± 17) x 10" 11 . 
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TABLE I. Event numbers of the e + e~ — > 7r + 7r~7r°(7) process (after background 
subtraction) and N bkg of background processes, integrated luminosity IL and detection 
efficiency e{s,E 1 — 0) (without 7-quantum radiation). S ra d is the radiative correction 
[Srad = £{s)/e(s,Ey = 0), £(s) is defined via the expression (19)]. 



^/s (MeV) 


T T / 1 — 1 \ 

IL (nb ) 




e(s,£ 7 = 


= o) 


N 3n 


Nbkg 


$rad 
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271.4±2.7 





,2544±0. 


.0082 


800± 34 


33±7 


0.905 


958 


249.0±2.5 





.2544±0 


.0082 


658± 29 


27±5 


0.918 


950 
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.2540±0. 


.0086 
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940 
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.0075 


1203± 41 


53±10 


0.937 


920 
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.2699±0. 


.0075 


1292± 42 


52± 8 
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58± 11 
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0.860 


660 


281.1±2.5 
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40± 11 
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120.1±1.2 





.3532±0. 
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0± 5 
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0.865 


600 


90.6±0.9 





3298±0. 


.0066 
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580 12.7±0.2 0.3561±0.0069 2± 4 2± 1 0.867 

560 11.2±0.2 0.3369±0.0068 -1± 1 1± 1 0.867 

540 12.1±0.2 0.3156±0.0067 -4± 2 4± 2 0.867 

520 7.2±0.2 0.2866±0.0065 0.861 

500 8.0±0.2 0.2278±0.0060 0± 1 1± 1 0.856 

480 13.4±0.2 0.2030±0.0058 0.852 

440 6.2±0.1 0.0183±0.0019 0.820 



TABLE II. The ratio of the cross sections, obtained by using different methods of 
background subtraction. is the number of the e + e~ — * 37r events obtained by fitting 
the two-photon invariant mass spectra. er(l) - the cross section measured in the approach 
when background was calculated according to Eq.(l), - the cross section measured by 
using the two-photon invariant mass spectra analysis. 



■sfs (MeV) 


N 3 n 




750 


1350±42 


0.995±0.034 


720 


700±32 


0.999±0.054 


690 


66±11 


0.956±0.209 


660 


17±10 


2.562±1.114 
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TABLE III. The e + e~ — > -k + iv~-k {) cross section. a mo d is the model uncertainty, atk g is 
the error due to background subtraction, a e ff © 07 l - the error due to uncertainty in the 
detection efficiency and integrated luminosity determination (3.4% &ty/s < 900 MeV and 
4.5% for^/i > 900 MeV), a sys = a e ff © cil © o- mod (s) © a bkg (s) is the total systematic 
error. 
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630.0 


< 0.34 (90% CL) 


600.0 


< 0.53 (90% CL) 


580.0 


< 1.36 (90% CL) 


560.0 


< 0.44 (90% CL) 


540.0 


< 1.21 (90% CL) 


520.0 


< 1.3 (90% CL) 


500.0 


< 0.96 (90% CL) 


480.0 


< 0.99 (90% CL) 


440.0 


< 24.7 (90% CL) 



TABLE IV. Fit results for the e + e~ — > 7r + 7r~7r° and lotv + -k~ cross sections. The column 
number N corresponds to the different models for the A pTr amplitude. Nf it is the number 
of fitted points. The first error is statistical, the second error shows the difference in the 
fit results due to various assumptions about the e + e~ — > 7r + 7r~7r° reaction dynamics and 
relative systematics between the SND and DM2 measurements. 



N 


1 


2 


3 


m u , MeV 


782.75±0.08 


782.72±0.08 


782.71±0.08 


T w , MeV 


8.60±0.04±0.01 


8.73±0.04±0.02 


8.63±0.04 


a(u> — ► 3tt), nb 


1609±7±1 


1624±10±5 


1610±7 ±3 


a((j) — > 3ir), nb 


645±6±4 


645±7±6 


658±8±7 


0^, degree 


161±2±4 


163±1±4 


187±4±6 


r , GeV- 1 






2.6± ( 1 ) i±0.2 


m u i, MeV 


1358±20±45 


1460±£§±70 


1410±30±60 


O, MeV 


500±|8±80 


1120±|88±200 


617±40±95 


a(uj' — > 37r), nb 


5.7±;j;|±0.5 


3.7±0.7±0.4 


5.0±0.2±0.4 


m u ii, MeV 


1808±Ig±20 


1760±50±40 


1750±20±6 


IV, MeV 


S07±i88±213 


540±?8g±50 


373±50±15 


a(u)" — > 37r), nb 


1.48±0.40±0.46 


2.4±0.7±0.9 


2.7±0.4±1.3 


a(uj" — > a;7r + 7r~), nb 


1.54±0.30±0.45 


1.8±0.4±0.5 


2.2±0.3±0.4 


a(p — > 37r), nb 




0.13±Ht±0.02 




0^p, degree 




-137±^±7 




xl/Nfit 
^i)/JV>« 

X^SND)/ N fit 
X3iv(DM2)/ N fit 


(80^120)/49 


(56^62)/49 


(45-^50)/49 


(21-j-46)/6 


(13.6-M6.4) /6 


(7.3-M1.3)/6 


(66-f-92)/67 


(49^56)/67 


(58^67)/67 


(22-j-37)/18 


(22-j-44)/18 


(27-=-42)/18 


xLn(DM2)/N 


(11-M5)/18 


10/18 


(23-7-26)/18 


Xlt/Nfit 


(182^260)/152 


(137-j-170)/152 


(156-180) /152 
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TABLE V. Fit results for the e + e~ -» tt+tt^tt and lott + tt' cross sections. The DM2 
data for e + e~ — > -k + tt~-k° were not used. The column number N corresponds to the 
different models for the A pw amplitude. Nfu is the number of fitted points. The first 
error is statistical, the second error shows the difference in the fit results due to various 
assumptions about the e + e _ — > 7r + 7r~7r° reaction dynamics. 



N 


1 


2 


3 


m„, MeV 


782.76±0.08 


782.72±0.08 


782.74±0.08 


r w , MeV 


8.63±0.04±0.0f 


8.7f±0.04±0.0f 


8.65±0.04 


a(u) — > 3ir), nb 


f608±7±f 


1618±10±2 


1612±7±2 


a((f) — * 3ir), nb 


653±6±4 


652±8±5 


665±10±8 


4>uj4„ degree 


166±3±3 


f65±3±4 


195±6±5 


r , GeV- 1 






2.3± l 1 ,;|±0.2 


m u i, MeV 


1273±|g±28 


f386±|jj±60 


1300±30±30 


IV, MeV 


405±||]±73 


827±lgg±186 


595±50±50 


a(u)' — ► 37r), nb 


6.9±0.4±0.7 


5.0±1.0±0.5 


5.6±0.3±0.5 


m u „, MeV 


1819±|8±32 


1773±sg±12 


1758±20±5 


I>, MeV 


679±ig±121 


505±i§!]±35 


345±50±10 


a(io" — > 37r), nb 


5.6±2.0±1.1 


5.7±1.7±0.6 


7.6±1.6±1.6 


cr(a>" — > o;7r + 7r _ ), nb 


1.2±0.3±0.2 


1.5±0.2±0.1 


1.7±0.2±0.1 


cr(p — > 37r), nb 




0.083±g;gi±0.009 




degree 




-134±il±8 




x|i)/JV/« 


(60^63)/49 


(52-^55)/49 


(40-^42)/49 


(10±16.6)/6 


(11.9±13)/6 


(4.6±6.3)/6 


X(SND)/Nfit 


(69^74)/67 


(51-=-52)/67 


(52^56)/67 


X3Tr(DM2)/Nfit 








Xwnn(DM2)/N 


(11±14)/18 


10/18 


23/18 


Xtot/Nfu 


(139±149)/134 


(112-^118)/ 134 


(115-M20) /134 



TABLE VI. Comparison of the ui — > 37r, it 7 and e + e decays branching ratios obtained 
by using the SND data with the world averages 1 . 





SND 


PDG-2002 


B(w -> e+e-) 
B(oj -> 3?r) 

B(0J -> 7T°7) 


(7.52 ±0.24) x 10~ b 
0.8965 ± 0.0051 
0.0865 ± 0.0045 


(6.95 ±0.15) x lO" 5 
0.8910 ±0.007 
0.087 ±0.004 
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FIG. 4. The # distribution of charged pions from the reaction e + e — > 7r + 7r 7r°. 
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FIG. 7. The angular distribution of photons. 
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FIG. 9. The photon energy distribution. 
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FIG. 10. The charged pions energy distribution 
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FIG. 12. The experimental xiir distribution in the energy regiony's > 870 MeV, fitted 
by a sum of distributions for the signal and background. The background contribution is 
shown by the filled histogram. 
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FIG. 13. The two-photon invariant mass distribution aty^ = 720 MeV. 
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FIG. 14. The detection efficiency e(E 1 ) dependence on the radiated photon energy 
for the e + e~ — + 7r + 7r~7r°(7) events at^/s ~ m w , obtained by simulation. 
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FIG. 17. The e + e~ — > 7r + 7r~7r° cross section. Dots are the SND data obtained in this 
work. Curves are results of fitting to the data in the model 2 (solid curve) and in the 
model 3 (dashed curve). 
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FIG. 18. The e + e~ — > 7r + 7r~7r° cross section. Dots are the SND data obtained in this 
work; the curve is the fit result. 
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FIG. 19. The e + e _ — > tt + -k~-k cross section. Dots are the SND data obtained in this 
work; the curve is the fit result. 
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FIG. 20. The e + e~ — > tt + tt~-k cross section. Dots are the SND data obtained in this 
work and in Ref. 27,2 . Curves are results of fitting to the data in the model 2 (solid curve) 
and in the model 3 (dashed curve). 
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FIG. 21. The e + e — > -k + tx it cross section. Dots are the SND data obtained in Ref. 27 ; 
the curve is the fit result. 
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FIG. 23. The e+e" -> wtt+tt" cross section. The results of the DM2 25 and CMD2 41 
experiments are shown. Curves are results of fitting to the DM2 data in the model 2 (solid 
curve) and in the model 3 (dashed curve) . 
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FIG. 24. The e+e~ -> tt+tt^tt cross section. The results of the SND 27 ' 2 and DM2 25 are 
shown. Curves are results of fitting to the data in the model 2. Dashed curve corresponds 
to the fit under assumption that a relative bias between the SND and DM2 data exists 
(DM2 data were scaled by a factor of C = 1.54). Solid curve is the result of the fitting to 
the SND data only. 
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FIG. 25. The ratio of the e + e — > 7r + 7r tt cross section obtained in different experi- 
ments to the fit curve. The shaded area shows the systematic error of the SND measure- 
ments. The SND (this work), DM1 17 , ND 8 ' 19 , CMD2 10 results are presented. 
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FIG. 26. The ratio of the e + e — > n + -K n cross section obtained in different exper- 
iments to the fit curve. The shaded area shows the systematic error of the SND mea- 
surements. The SND (this work and Ref. 27 ' 2 ) DM1 17 , ND 8 ' 19 , CMD2 15 ' 16 results are 
presented. 
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FIG. 27. The ratio of the e + e — > 7r + 7r 7r cross section obtained in different experi- 
ments to the fit curve. The shaded area shows the systematic error of the SND measure- 
ments. The SND 27 , DM1 17 , ND 8,19 , CMD2 15 ' 16 results are presented. 
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FIG. 28. The ratio of the e e — ► 7r ir 7r cross section obtained in different experi- 
ments to the fit curve. The shaded area shows the systematic error of the SND measure- 
ments. The SND 27 ' 2 , ND 8 ' 19 , CMD2 16 results are presented. 
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FIG. 29. The w meson mass measured in this work 

Ref. 10 ' 42 ' 43 ' 6 ' 44 ' 17 . The shaded area shows the world average value 1 . 
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FIG. 30. The u meson width r u measured in this work (SND-03) and in 
Ref. 10,42 ' 7 ' 9,6 ' 17,5 . The shaded area shows the world average value 1 . 
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FIG. 31. The value of a{tu -> 3vr) measured in this work (SND-03) and in R e f. 10 . 7 . 9 . 6 > 17 . 5 . 
The shaded area shows the world average value 1 . 
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FIG. 32. The ratio of the partial widths lo —* ty°j and uj — > 37r, obtained in this work 
(SND-03) and in Ref. 45 ' 7 ' 44 ' 5 ' 46 . The shaded area shows the world average value 1 . 
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FIG. 33. The value of a(<f) — > 37r) obtained in this work (SND-03) and in 
Ref. 16 ' 15 ' 8 ' 18 ' 17 ' 14,12,11 . The shaded area shows the world average value according to the 
year 2000 PDG table 48 . 
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FIG. 34. The phase 0^ obtained in this work. The star and the dot indicate the phase 
values obtained from the fit in the second and third models respectively. The shaded areas 
show the expected energy behavior of the 4>u><t> phase 39 for various values of the range 
parameter ro. 
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